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We study the fault tolerance of quantum computation by adiabatic evolution, a quantum al- 
gorithm for solving various combinatorial search problems. We describe an inherent robustness of 
adiabatic computation against two kinds of errors, unitary control errors and decoherence, and we 
study this robustness using numerical simulations of the algorithm. 
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I. INTRODUCTION 

The method of quantum computation by adiabatic 
evolution has been proposed as a general way of solv- 
ing combinatorial search problems on a quantum com- 
puter jjj. Whereas a conventional quantum algorithm 
is implemented as a sequence of discrete unitary trans- 
formations that form a quantum circuit involving many 
energy levels of the computer, the adiabatic algorithm 
works by keeping the state of the quantum computer 
close to the instantaneous ground state of a Hamiltonian 
that varies continuously in time. Therefore, an imperfect 
quantum computer implementing a conventional quan- 
tum algorithm might experience different sorts of errors 
than an imperfect adiabatic quantum computer. In fact, 
we claim that an adiabatic quantum computer has an in- 
herent robustness against errors that might enhance the 
usefulness of the adiabatic approach. 

The adiabatic algorithm works by applying a time- 
dependent Hamiltonian that interpolates smoothly from 
an initial Hamiltonian whose ground state is easily pre- 
pared to a final Hamiltonian whose ground state encodes 
the solution to the problem. If the Hamiltonian varies 
sufficiently slowly, then the quantum adiabatic theorem 
guarantees that the final state of the quantum computer 
will be close to the ground state of the final Hamiltonian, 
so a measurement of the final state will yield a solution 
of the problem with high probability. This method will 
surely succeed if the Hamiltonian changes slowly. But 
how slow is slow enough? 

Unfortunately, this question has proved difficult to an- 
alyze in general. Some numerical evidence suggests the 
possibility that the adiabatic algorithm might efficiently 
solve computationally interesting instances of hard com- 
binatorial search problems, outperforming classical meth- 
ods Whether the adiabatic algorithm provides a 
definite speedup over classical methods remains an in- 
teresting open question. As we will discuss in Sec. g 
the time required by the algorithm for a particular in- 



stance can be related to the minimum gap A between 
the instantaneous ground state and the rest of the spec- 
trum. Roughly speaking, the required time goes like 
A~ 2 . Thus, if A~ 2 increases only polynomially with the 
size of the problem, then so does the time required to 
run the algorithm. However, determining A has not been 
possible in general. 

Our objective in this paper is not to explore the com- 
putational power of the adiabatic model, but rather to 
investigate its intrinsic fault tolerance. Since quantum 
computers are far more susceptible to making errors than 
classical digital computers, fault tolerant protocols will 
be necessary for the operation of large-scale quantum 
computers. General procedures have been developed that 
allow any quantum algorithm to be implemented fault 
tolerantly on a universal quantum computer || , but these 
involve a substantial computational overhead. Therefore, 
it would be highly advantageous to weave fault tolerance 
into the design of our quantum hardware. 

We therefore will regard adiabatic quantum computa- 
tion not as a convenient language for describing a class of 
quantum circuits, but as a proposed physical implemen- 
tation of quantum information processing. We do not 
cast the algorithm into the conventional quantum com- 
puting paradigm by approximating it as a sequence of 
discrete unitary transformations acting on a few qubits 
at a time. Instead, suppose we can design a physical de- 
vice that implements the required time-dependent Hamil- 
tonian with reasonable accuracy. We then imagine im- 
plementing the algorithm by slowly changing the param- 
eters that control the physical Hamiltonian. How well 
does such a quantum computer resist decoherence, and 
how well does it perform if the algorithm is imperfectly 
implemented? 

Regarding resistance to decoherence, we can make a 
few simple observations. The phase of the ground state 
has no effect on the efficacy of the algorithm, and there- 
fore dephasing in the energy eigenstate basis is presum- 
ably harmless. Only the interactions with the environ- 
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ment that induce transitions between eigenstates of the 
Hamiltonian might cause trouble. In principle, these can 
be well controlled by running the algorithm at a tem- 
perature that is small compared to the minimum gap A. 
(We use units in which Boltzmann's constant hp, = 1 ; 
so that temperature has units of energy.) If A decreases 
slowly as the size of the problem increases, then the re- 
sources required to run at a sufficiently low temperature 
may be reasonable. Since the adiabatic method is only 
efficient if A is not too small, we conclude that whenever 
the method works on a perfectly functioning quantum 
computer, it is robust against decoherence. 

In addition to environmental decoherence, we must 
also consider the consequences of imperfect implemen- 
tation. Our chosen algorithm may call for the time- 
dependent Hamiltonian H(t), but when we run the al- 
gorithm, the actual Hamiltonian will be H(t) + Kit), 
where K (t) is an "error." An interesting feature of adia- 
batic quantum computation is that K(t) need not remain 
small during the evolution in order for the algorithm to 
work effectively. A reasonably large excursion away from 
the intended Hamiltonian is acceptable, as long as K it) 
is slowly varying and has initial and final values that are 
not too large. A very rapidly fluctuating K(t) may also 
be acceptable, if the characteristic frequency of the fluc- 
tuations is large compared to the energy scale of H(t). 

In this paper, we use numerical simulations to inves- 
tigate the sensitivity of an adiabatic computer to deco- 
hering transitions and to a certain class of unitary per- 
turbations induced by a Hamiltonian Kit). The results 
are consistent with the idea that the algorithm remains 
robust as long as the temperature of the environment 
is not too high and K (t) varies either sufficiently slowly 
or sufficiently rapidly. Thus, the adiabatic model illus- 
trates the principle that when the characteristics of the 
noise are reasonably well understood, it may be possible 
to design suitable quantum hardware that effectively re- 
sists the noise. However, note that some of the effects 
of decoherence and unitary control error may not be sig- 
nificant for the small problems we are able to study — 
especially in the case of decoherence, where the time re- 
quired by the simulation restricts us to systems with only 
four qubits — and hence our data may not be indicative 
of the performance of the algorithm working on larger 
inputs. 

A technique closely related to adiabatic computation 
was described by Kadowaki and Nishimori Q and has 
been tested experimentally (in conjunction with a cooling 
procedure) by Brooke et al. 0. In a different guise, the 
principles that make quantum adiabatic evolution robust 
also underlie the proposal by Kitaev || to employ non- 
abelian anyons for fault-tolerant quantum computation. 
The fact that adiabatic evolution incorporates a kind of 
intrinsic fault tolerance has also been noted in [f^bl . 

In Sec. g we review the adiabatic model of quantum 
computation, and in Sec. Ill we describe the specific com- 
binatorial search problem (three-bit exact cover) that we 
use in our simulations. Sees. IV and |y| report our nu- 



merical results on decoherence and unitary control error, 
and Sec. VI summarizes our conclusions. 



II. ADIABATIC QUANTUM COMPUTATION 

We briefly review the adiabatic model of quantum com- 
putation introduced in Let h(z) be a function of n 
bits z — (z\, z%, Z3, . . . , z n ), and consider the computa- 
tional problem of finding a value of z that minimizes h(z) . 
We will typically be interested in the case were this value 
of z is unique. We may associate with this function the 
Hermitian operator 
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h p =J2 M*)I*X*I> 



(1) 



z=0 



so that the computational basis state \z) is an eigenstate 
of Hp with eigenvalue h(z). Then the problem is to de- 
termine which state \z) is the ground state (eigenstate 
with lowest eigenvalue) of Hp. We refer to Hp as the 
problem Hamiltonian. 

The strategy for finding the ground state of Hp is to 
prepare the ground state of some other beginning Hamil- 
tonian Hb and slowly interpolate to Hp . A simple choice 
for the interpolation is given by the one-parameter family 
of Hamiltonians 



H{s) = (1 - s)H B + sH P 



(2) 



that interpolates between Hp and Hp as s varies from 
to 1. We prepare the ground state of Hp at time t = 0, 
and then the state evolves from t = to t = T according 
to the Schrodinger equation, 



where the Hamiltonian is 



H[t) = Hit/T) . 



(3) 



(4) 



At time T (the run time of the algorithm), we measure 
the state in the computational basis. If we let \<p) denote 
the (unique) ground state of Hp for a given instance of 
the problem, then the success probability of the algorithm 
for this instance is 



Prob(T)^|(y#(T)}| 2 



(5) 



Does the algorithm work? According to the quantum 
adiabatic theorem |l5|,[l6|], if there is a nonzero gap be- 
tween the ground state and the first excited state of H(s) 
for all s € [0, 1], then Prob(T) approaches 1 in the limit 
T — > 00. Furthermore, level crossings are non-generic 
in the absence of symmetries, so a non- vanishing gap is 
expected if Hp does not commute with Hp. Thus, the 
success probability Prob(T) of the algorithm will be high 
if the evolution time T is large enough. The question is: 
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how large a T is large enough so that Prob(T) is larger 
than some fixed constant? 

We can reformulate this question in terms of 



A = min (Ei(s) - E (s)) 
se[o,i] 



and 



max 

sG[0,l] 



, AH, v 

;mi^io, s) 



(6) 



(7) 



where Eq(s) is the lowest eigenvalue of H(s), Ei(s) is the 
second lowest eigenvalue, and 0, s) , |l,s) are the corre- 
sponding eigenstates. By calculating the transition prob- 
ability to lowest order in the adiabatic expansion [jL6j, 
one finds that the probability of a transition from ground 
state to first excited state is small provided that the run 
time T satisfies 



(8) 



If the spectrum consists of only two levels, then this con- 
dition is sufficient to ensure that the system remains in 
the ground state with high probability. In general, the 
required run time T will be bounded by a polynomial in 
n so long as A and £ are polynomially bounded. For the 
problems we are interested in, £ is polynomially bounded, 
so we only have to consider the behavior of A. 

By rescaling the time, we can think of the evolution 
as taking place in the unit time interval between s = 
and s = 1, but in that case the energy eigenvalues are 
rescaled by the factor T. Roughly speaking, we can 
think of dH(s)/ds as a perturbation that couples the 
levels of the instantaneous Hamiltonian H(s), and has 
the potential to drive a transition from \0,s) to \l,s). 
But if T is large, the effects of this perturbation are 
washed out by the rapid oscillations of the relative phase 
expH'T JSds'iE^s') - E (s'))}. 

Note that the Hamiltonian can be regarded as reason- 
able only if it is "local," that is, if it can be expressed 
as a sum of terms, where each term acts on a bounded 
number of qubits (a number that does not grow with n). 
Indeed, in this case, the Hamiltonian evolution can be 
accurately and efficiently simulated by a universal quan- 
tum computer p7| . Many combinatorial search problems 
(e.g., 3SAT) can be formulated as a search for a minimum 
of a function that is local in this sense. Along with a lo- 
cal choice of Hb, this results in a full H(t) that is also 
local. 

A direct physical implementation of the continuously 
varying H (t) would presumably be possible only under a 
somewhat stronger locality condition. We might require 
that each qubit is coupled to only a few other qubits, 
or perhaps that the qubits can be physically arranged in 
such a way that the interactions are spatially local. For- 
tunately, there are interesting computational problems 
that have such forms, such as 3SAT restricted to having 



each bit involved in only three clauses or the problem of 
finding the ground state of a spin glass on a cubic lat- 
tice [ H . However, for the purposes of our simulation, 
we will only consider small instances, and since we do 
not have a specific physical implementation in mind, we 
will not concern ourselves with the spatial arrangement 
of the qubits. 



III. AN EXAMPLE: THE EXACT COVER 
PROBLEM 

For definiteness, we study the robustness of the adi- 
abatic algorithm via its performance on the problem 
known as "three-bit exact cover" (EC3). An n-bit in- 
stance of EC3 consists of a set of clauses, each of which 
specifies three of the n bits. A clause is said to be satis- 
fied if and only if exactly one of its bits has the value 1. 
The problem is to determine if any of the 2 n assignments 
of the n bits satisfies all of the clauses. 

For this problem, the function h(z) is a sum 



H z ) = ^h c {z ic ,z jc ,z kc ) 



(9) 



of three-bit clauses, where 



he { Z ic J z jc > z kc ) 

, (zi c , zj c , Zk c ) satisfies clause C , » , 
Zi c , Zj c , Zk c ) violates clause C . 

The value of the function h(z) is the number of clauses 
that are violated; in particular, h(z) — if and only if z 
is an assignment that satisfies all the clauses. 

To solve EC3 by the adiabatic algorithm, a sensible 
choice for the beginning Hamiltonian is 



where 



Hb — E Hb,c 



H B ,a = l(l-^ ) ) + l(l-4 io) 



(11) 



4(1-^)' 
which has the ground state 



|V(0)> 



2™/ 2 



(12) 



(13) 



The resulting H(t) is local in the sense that it is a sum 
of terms, each of which acts on only a few qubits. A 
stronger kind of locality can be imposed by restricting 
the instances so that each bit is involved in at most a 
fixed number of clauses. The computational complexity 
of the problem is unchanged by this restriction. 
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Numerical studies of the adiabatic algorithm applied 
to this problem were reported in J|,|| . Instances of EC3 
with n bits were generated by adding random clauses 
until there was a unique satisfying assignment, giving 
a distribution of instances that one might expect to be 
computationally difficult to solve. The results for a small 
number of bits (n < 20) were consistent with the possi- 
bility that the adiabatic algorithm requires a time that 
grows only as a polynomial in n for typical instances 
drawn from this distribution. If this is the case, then the 
gap A does not shrink exponentially. Although the typ- 
ical spacing between levels must be exponentially small, 
since there are an exponential number of levels in a poly- 
nomial range of energies, it is possible that the gap at the 
bottom is larger. For example, Fig. ^ shows the spectrum 
of a randomly generated seven-bit instance of EC3. The 
gap at the bottom of the spectrum is reasonably large 
compared to the typical spacing. This feature is not spe- 
cific to this one instance, but is characteristic of randomly 
generated instances, at least for n < 10, beyond which 
the repeated matrix diagonalization required to create a 
picture of the spectrum becomes computationally costly. 
A large gap makes an instance readily solvable by the 
adiabatic algorithm, and also provides robustness against 
thermal transitions out of the ground state. 




FIG. 1. Spectrum of a randomly generated n = 7 bit in- 
stance of EC3 with a unique satisfying assignment. Note that 
the energy gap between the ground state and the first excited 
state is significantly larger than all other gaps. An expanded 
view would show that there are no level crossings anywhere 
in the spectrum (except for the degeneracies at s = and 
s = l). 



IV. DECOHERENCE 

Perhaps the most significant impediment to building a 
large-scale quantum computer is the problem of decoher- 
ence. No quantum device can be perfectly isolated from 



its environment, and interactions between a device and 
its environment will inevitably introduce noise. Fortu- 
nately, such effects can be countered using fault-tolerant 
protocols, but as we have already mentioned, these pro- 
tocols can be costly. Therefore, we would like to consider 
quantum systems with inherent resistance to decohering 
effects. If the ground state of our adiabatic quantum 
computer is separated from the excited states by a sizable 
energy gap, then we expect it to exhibit such robustness. 
Here, we consider how the adiabatic algorithm for EC3 
is affected by decoherence. 

First, we briefly review the master equation formalism 
for describing the decohering effects of an environment 
on a quantum system. Suppose that our quantum com- 
puter is a collection of spin-i particles interacting with 
each other according to the Hamiltonian H$ and weakly 
coupled to a large bath of photons. The total Hamilto- 
nian of the quantum computer and its environment is 



H = Hc 



XV. 



(14) 



where He is the Hamiltonian of its environment, V is 
an interaction that couples the quantum computer and 
the photon bath, and A is a coupling constant. We may 
describe the state of the quantum computer alone by the 
density matrix p found by tracing over the environmen- 
tal degrees of freedom. In general, the time evolution of 
p is complicated, but under reasonable assumptions, we 
can approximate its evolution using a Markovian master 
equation. 

One way of deriving such a master equation is to con- 
sider the weak coupling limit, in which A <C 1 If the 
environment is very large and only weakly coupled to the 
quantum computer, it will be essentially unchanged by 
the interaction. Furthermore, in this limit, we can expect 
the evolution of the quantum computer to be Markovian, 
or local in time, if we filter out high-frequency fluctua- 
tions by some coarse-graining procedure. Assuming that 
the combined state of the quantum computer and its en- 
vironment begins in a product state p(0) ® pe, Davies 
derives the master equation 



dp 
dt 



= -i{H s ,p]+\ 2 K*p, 



(15) 



where 

Kp = - &xTt e [U{-x)VU{x),[V,p® p E }} (16) 
Jo 

K^p= lim - / dyU(-y){K[U(y)pU(-y)]}U(y) (17) 



with 



U(x) 



-ix(H s +H E ) 



(18) 



where we have (temporarily) assumed that H$ is time- 
independent. Although the tj operation defined by ( |l7j ) 
does not appear in some formulations of the Markovian 
master equation, it appears to be essential for the equa- 
tion to properly describe the weak coupling limit [2(1, 
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and in particular, for it to capture the physics of relax- 
ation to thermal equilibrium. The master equation (|l5| ) 
has the property that if the environment is in thermal 
equilibrium at a given temperature, then the decoher- 
ing transitions drive the quantum computer towards the 
Gibbs state of Hs at that temperature. While not an 
exact description of the dynamics, ( plf ) should provide a 
reasonable caricature of a quantum computer in a ther- 
mal environment. 

Note that (|ll]) is derived assuming a time-independent 
Hamiltonian Hs; with a time- varying Hs(t), we should 
expect the generator of time evolution at any particu- 
lar time to depend on the Hamiltonian at all previous 
times ^J. However, if Hs(t) is slowly varying, then it is 
a good approximation to imagine that the generator at 
any particular time depends only on Hs at that time [p2[ . 
In particular, since we are interested in nearly adiabatic 
evolution, Hs{t) varies slowly, so (|l^) remains a good ap- 
proximation, where at any given time t we compute 
using only Hs(t). Note that with Hs(t) time-dependent, 
U(x) defined by ( |l8|) is not the time evolution operator; 
it depends on the time t only implicitly through Hs(t). 

For a system of spins coupled to photons, we choose 
the interaction 



V = 



du>[g(u)a u a¥ +g*(u)aia^}, (19) 



where is a sum over the spins, are raising and 
lowering operators for the ith spin, a u is the annihila- 
tion operator for the photon mode with frequency u>, and 
\g{ui) is the product of the coupling strength and spectral 
density for that mode. Note that if the coupling strength 
is frequency-dependent, we can absorb that dependence 
into g(w), leaving A as a frequency-independent parame- 
ter. With this specific choice for V, we can perform the 
integrals and trace in (|l5|-]l8|) . If we assume that all spac- 
ings between eigenvalues of Hs are distinct, the resulting 
expression simplifies considerably, and we find 



dp 
dt 



-i[H s , P ] 

-^[Nb*\9b*\ 2 {aW®\b)(b\a%) 

i,a,b 

+ (N ab + l)\g ab \ 2 {b\a ( l ) \a}{a\cTf\b}} 
{(| )(a|p) + (p|a)<a|)-2|6)( Ho>(6|} > (20) 



where the states \a) are the time-dependent instanta- 
neous eigenstates of Hs with energy eigenvalues w a , 



N ba 



1 



(21) 



exp[/3(w 6 - LOa)] - 1 

is the Bose-Einstein distribution at temperature 1/(3, and 



9ba 



Xg(cob - 0J a 
, 



bJb > UJa 
LO b <U} a . 



We simulated the effect of thermal noise by numerically 
integrating the master equation with a Hamiltonian 
Hs given by (|]) and with the initial pure state density 
matrix p(0) — \ip(0))(ip(0)\ given by (|15|). For simplicity, 
we chose g(u>) = 1 for to > and zero otherwise. Al- 
though we would expect that g(to) — » as to — ► oo, for 
the small systems we are able to simulate it should be a 
reasonable approximation to treat g{to) as constant and 
tune the overall coupling strength using A 2 . 

How should we expect the success probability 
(tp\p(T)\(p) , where \tp) is the ground state of Hp, to de- 
pend on the run time T and the temperature? If the run 
time T is sufficiently long, then regardless of its initial 
state the quantum computer will come to thermal equi- 
librium; at the time of the final readout it will be close 
to the Gibbs state 



e -0H P 
lim p(T) = — wtz- 



PP 



(23) 



(22) 



of the problem Hamiltonian Hp, and the success proba- 
bility will be approximately (tp\pp\<p). This probability 
may be appreciable if the temperature is small compared 
to the gap between the ground state and first excited 
state of Hp. Thus one way to find the ground state of 
Hp is to prepare the computer in any initial state, put 
it in a cold environment, wait a long time, and measure. 
However, this thermal relaxation method is not an effi- 
cient way to solve hard optimization problems. Although 
it may work well on some instances of a given problem, 
this method will not work in cases where the computer 
can get stuck in local minima from which downward tran- 
sitions are unlikely. In such cases, the time for equilibra- 
tion is expected to be exponentially large in n. 

Consider an instance with a long equilibration time 
so that cooling alone is not an efficient way to find the 
ground state of Hp. It is possible that the minimum gap 
A associated with the quantum algorithm is not small, 
and the idealized quantum computer, running without 
decohering effects, would find the ground state of Hp in 
a short time. In this situation, if we include the cou- 
pling of the system to the environment and we run at a 
temperature much below A, then thermal transitions are 
never likely, and the adiabatic algorithm should perform 
nearly as well as in the absence of decoherence. But if the 
temperature is comparable to A, then the performance 
may be significantly degraded. 

On the other hand, consider an instance for which the 
equilibration time is short, so that cooling alone is a good 
algorithm. Furthermore, suppose that the adiabatic al- 
gorithm would find the ground state of Hp in a short 
time in the absence of decohering effects. In this case, 
the combined effects of cooling and adiabatic evolution 
will surely find the ground state of Hp in a short time. 
But note that A alone does not control the success of 
the algorithm. Even if H(t) changes too quickly for the 
evolution to be truly adiabatic so that a transition occurs 
where the gap is smallest, the system may be cooled back 
into its ground state at a later time. 
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FIG. 2. The success probability as a function of run time T for two instances of EC3 with n = 4 bits. The instance on the 
left has a gap of Ai ~ 0.301 and the instance on the right has a gap of A2 « 0.425. The dotted line shows the behavior of 
the algorithm with no decoherence, i.e., A 2 = 0. Note that in the figure on the right, the dotted curve is partially obscured 
but can be seen slightly above the topmost solid curve. The solid lines show the behavior of the algorithm in the presence of 
decoherence with A 2 = 0.1 for five different temperatures. The triangles at the far right show the thermal success probabilities 
(<p\pp\<p) at each of these temperatures. From top to bottom, the temperatures are 1/10, 1/2, 1, 2, and 10. 



Typical results of the simulation are shown in Fig. ^ 
for two 71 = 4 bit instances of EC3 with unique satis- 
fying assignments. These two instances have minimum 
gaps of Ai « 0.301 and A2 « 0.425. For each instance, 
we plot the success probability as a function of the run 
time T. With A 2 = 0.1, we consider five temperatures: 
1/10, 1/2, 1, 2, and 10. We also present the data with 
no decoherence (A 2 = 0) for comparison. 

Unfortunately, the time required to integrate j20| ) 
grows very rapidly with n. Whereas a state vector con- 
tains 2™ entries, the density matrix contains 4™ entries; 
and in addition, calculating dp/dt at each timestep re- 
quires evaluating a double sum over 2™ energy eigen- 
states. For this reason, we were only able to consider 
instances with n < 4. 

The results are consistent with our general expecta- 
tions. In the absence of decoherence, the success proba- 
bility becomes appreciable for sufficiently long run times. 
This probability rises faster for the problem with a larger 
gap. When we add decoherence at high temperature, the 
success probability never becomes very large (note the 
lowest curves in Fig. |J) . As the temperature is decreased 
to a value of order one, the presence of decoherence has a 
less significant effect on the success probability. In fact, 
for sufficiently low temperatures, the success probabil- 
ity can actually be higher in the presence of decoher- 
ence than when there is no decoherence. This is because 
the primary effect of decoherence at low temperature is 
to drive transitions towards the ground state, improving 
performance. 

However, these results do not illustrate a definitive con- 
nection between the minimum gap A and the tempera- 
ture above which the algorithm no longer works. These 
simple n = 4 bit instances fall into the second category 



discussed above: the equilibration time is short, so cool- 
ing alone is a good algorithm. In other words, no sharp 
distinction can be drawn between the run time required 
for the adiabatic algorithm to perform well in the absence 
of decoherence and the run time required for equilibra- 
tion. Accordingly, the dependence of the success prob- 
ability on temperature and run time is similar for the 
two instances shown in Fig. |^, even though the minimum 
gaps for these instances are somewhat different. 



V. UNITARY CONTROL ERROR 

We now consider how the performance of the adiabatic 
algorithm for EC3 is affected by adding three different 
kinds of perturbations to the Hamiltonian. Each pertur- 
bation we consider is a sum of single-qubit terms, where 
each term can be interpreted as a magnetic field point- 
ing in a random direction. To simplify our analysis, we 
assume that the magnitude of the magnetic field is the 
same for all qubits, but its direction varies randomly from 
qubit to qubit. The perturbations we consider are 



Kt{a) 



K*(a) 



m* ■ a 



i=l 

n 

C2 sin("7rs) rhi ■ c?W , 
i=l 

1 " 

- sin(C 3 7rs) ^ ™i ' ^ W 1 



(24) 
(25) 
(26) 



i=l 



which are added to (E 
tonian according to (J 



and give a time-dependent Hamil- 
1. Each rhi is a randomly generated 



G 



real three-component vector with unit length, C\ and C2 
are real numbers, and C3 is a nonnegative integer. 

The adiabatic algorithm was simulated by numeri- 
cally solving the time-dependent Schrodinger equation 
with initial state \">p(0)) given by Jl3| ) and Hamiltonian 
H(t/T)+Kj(t/T) for a given j e {1,2,3}. As in §-|, 
we used a fifth-order Runge-Kutta method with variable 
step-size, and checked the accuracy by verifying that the 
norm of the state was maintained to one part in a thou- 
sand. For a specified value of n, we randomly generated 
an instance of EC3 with a unique satisfying assignment. 
Then we randomly generated several different values of 
the magnetic field directions {mi}. For each instance of 
the problem and the magnetic field, the run time was 
chosen so that the success probability without the per- 
turbation was reasonably high. With this run time fixed, 
we then determined the success probability for varying 
values of the relevant Cj. 

First, we consider the perturbation K\. Since it turns 



on at a constant rate, this perturbation can be thought 
of as an error in Hp. Note that with C\ 7^ 0, the final 
Hamiltonian is not simply Hp, so the algorithm will not 
work exactly even in the adiabatic limit T — > 00. This 
perturbation is potentially dangerous because of the way 
its effect scales with the number of bits n. Indeed, con- 
sider the case where Hp can be separated into a sum of 
Hamiltonians acting separately on each qubit. If adding 
K\ reduces the overlap of the ground state \<p) of Hp with 
the perturbed ground state \ip') by some fixed value e for 
each of the n qubits, then the total overlap is (1 — e) n , 
which is exponentially small in the number of bits. Thus 
the algorithm clearly fails in this factorized case. In gen- 
eral, if the magnitude of K\ is independent of n, then we 
expect the algorithm to fail. However, if the magnitude 
of K\ falls as 1/n or faster, then the shift of the ground 
state may be small enough (as it would be in the factor- 
ized case) that the algorithm is not significantly affected. 
Note that for any n there is some value of C\ that is small 
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FIG. 3. (Top) The success probability of the adiabatic algorithm for two randomly generated instances of EC3 with n = 7 
bits (left) and n = 10 bits (right) under the perturbation Ki defined by (|2^) for four different sets of magnetic field directions. 
For each n, the run time is the same for each random perturbation. (Bottom) The corresponding overlaps | (ip\p') | 2 of the 
ground state \tp) of Hp with the perturbed ground state \<p') at s = 1. 
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enough that the disadvantage of reduced overlap with the 
ground state of Hp may be overcome if the perturbation 
happens to increase the minimum gap A. For this reason, 
we expect to sometimes see an increase in success prob- 
ability for small C\ that goes away as C\ is increased. 

The effect of the perturbation K\ is shown in Fig. || 
for n = 7 and n = 10 bit instances of EC3, with four 
different randomly generated sets of magnetic field direc- 
tions for each instance. The run time is chosen such that 
for C\ = 0, the success probability is around 1/2. The 
top plots show that for small C\, the success probabil- 
ity is not strongly suppressed; in fact, in some cases it 
is significantly enhanced. For large enough Ci, the suc- 
cess probability is heavily suppressed. The bottom plots 
show the overlap |(c/?|<//)| 2 between the ground state of 
Hp and the actual ground state in the presence of the 
perturbation. As we expect, the suppression of the suc- 
cess probability is correlated with the amount of overlap. 
We also studied a similar perturbation in which s is re- 
placed by 1 — s, which can be thought of as an error in 
Hp. Unsurprisingly, the results were qualitatively simi- 



lar. 

Next, we consider the low frequency perturbation Ki. 
The period of oscillation is chosen such that the perturba- 
tion vanishes at t — and t = T, so the perturbation does 
not affect the algorithm in the adiabatic limit. Since the 
success probability is quite sensitive to the value of the 
minimum gap A, and it is not a priori obvious whether 
a perturbation will increase or decrease A, we can guess 
that turning on a nonzero value of C2 can either increase 
the success probability or decrease it. 

Fig. [| shows the effect of the perturbation K2, using 
the same instances, magnetic field directions, and run 
times as in Fig. ||. The top plots show the success prob- 
ability as a function of C2. As in the case of Ky, some 
perturbations can raise the success probability and some 
suppress it. Perhaps unsurprisingly, a particular set of 
magnetic field directions that can raise the success prob- 
ability under K\ is also likely to help when K2 is applied. 
But unlike K±, K2 can improve the success probability 
even with C2 — 2, where the size of the perturbation 
is comparable to the size of the unperturbed Hamilto- 
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nian. The bottom plots show the minimum gap A when 
the perturbation is added. Note that there is a strong 
correlation between the success probability and A. 

For both perturbations K\ and Ki , similar results have 
been observed (with fewer data points) for instances with 
as many as n = 14 bits. Figs. [| and ^ present typical 
data. For example, for a given instance, typically one or 
two out of four sets of randomly chosen magnetic field 
directions led to an improvement in the success proba- 
bility for some values of C\ and C2, compared to the 
unperturbed case. 

Finally, we consider the perturbation in which the 
magnitude of the oscillating component is fixed, but we 
may vary its frequency by varying C3. As for K2, the 
frequency is chosen so that the perturbation vanishes at 
t = and t = T. We expect that for C3 of order one, 
the perturbation will be likely to excite a transition, and 
that the success probability will be small. But since both 



Hb and Hp have a maximum eigenvalue of order n, we 
can anticipate that for 

nT 

C 3 »— , (27) 

7T 

the perturbation will be far from any resonance. Then 
the probability that the perturbation drives a transition 
will be low, and the success probability should be com- 
parable to the case where the perturbation vanishes. 

Some representative plots of the dependence of the suc- 
cess probability on C3 are shown in Fig. ||. Each plot 
corresponds to a particular randomly generated instance 
of EC3 (with either n = 8 bits or n = 10 bits) and a 
randomly generated set of magnetic field directions. In 
the top row of plots, the run time is chosen so that the 
success probability is around 1/8 with the perturbation 
absent (i.e., C3 = 0). In the bottom row, the run time is 



n = 8 n = 10 




C 3 C 3 

FIG. 5. The success probability as a function of the frequency C3 of the perturbation K3 defined in The data in each 

plot were obtained for a randomly generated instance of EC3 with randomly generated magnetic field directions. The data in 
the left column are for two instances with n = 8 bits, and the data in the right column are for two instances with n = 10 bits. 
For the top row, the run time is chosen so that the success probability is around 1/8 for C3 = 0, and for the bottom row, the 
run time is twice as long. The leftmost points in each plot correspond to C3 = 0, so the perturbation is absent for all t. C3 
takes integer values, so the lines are included only to guide the eye. 
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doubled. All of the data exhibit the expected qualita- 
tive trend. The leftmost point corresponds to C3 = 0. 
For the smallest values of C3 > 0, the success probabil- 
ity may not be too badly damaged; for somewhat larger 
values of C3 it is heavily suppressed; and for sufficiently 
large C3 it recovers to a value near the success probabil- 
ity in the absence of the perturbation. The value of hT/tt 
is around 19 and 39 for the upper and lower n = 8 plots 
and is around 38 and 76 for the upper and lower n = 10 
plots, so the estimate (|27]) turns out to be reasonable. 

Another conspicuous feature of the plots in Fig. || is 
that the success probability tends to oscillate between 
even and odd values of C3, though whether even or odd 
values are favored varies from case to case. This occurs 
because the perturbation's time average vanishes for C3 
even, so that its integrated effect is weaker than for C3 
odd. Since a small perturbation might either help or hurt, 
the success probability is slightly enhanced for odd C3 in 
some cases, and is slightly suppressed in other cases. 
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VI. CONCLUSIONS 



We have conducted numerical simulations to investi- 
gate the fault tolerance of adiabatic quantum computa- 
tion, and our results are consistent with the claim that 
this algorithm is robust against decoherence and certain 
kinds of random unitary perturbations. Thus, if a phys- 
ical system could be engineered with interactions rea- 
sonably well described by a Hamiltonian that smoothly 
interpolates from an initial Hp to a final Hp correspond- 
ing to an interesting combinatorial search problem, and if 
the gap remains large throughout the interpolation, that 
system might be a powerful computational device. 

Although we have viewed unitary perturbations as 
noise, the fact that they sometimes raise the success 
probability suggests a possible way to speed up the adi- 
abatic algorithm. The algorithm finds the ground state 
of Hp by starting the system in the ground state of Hb ■ 
The quantum state evolves as the system Hamiltonian 
smoothly interpolates from Hp to Hp. However, there 
are many possible choices for Hp and many smooth paths 
from a given Hp to Hp. The choices (^Tj) and (0) are 
convenient but arbitrary, so choosing an alternate route 
to Hp might speed up the algorithm. An example of 
this is seen in pj|, where it is shown that optimizing 
the time-dependent coefficients of Hp and Hp allows the 
adiabatic algorithm to achieve a square root speedup for 
an unordered search problem. More generally, the in- 
terpolating Hamiltonian might involve terms which have 
nothing to do with Hp or Hp, but which increase A and 
therefore improve performance. For example, the pertur- 
bation K2 sometimes increases the success probability, as 
seen in Fig. |^. Rather than being thought of as a source 
of error, such a perturbation could be applied intention- 
ally and might sometimes enhance the effectiveness of the 
adiabatic algorithm. 
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